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Jpt ' Abstract 

CN ; 

The XXZ spin chain with a boundary magnetic field h is considered, using 
the vertex operator approach to diagonalize the Hamihonian. We find explicit 
bosonic formulas for the two vacuum vectors with zero particle content. There 
tH- ' are three distinct regions when /i > 0, in which the structure of the vacuum 

CTS . states is different. Excited states are given by the action of vertex operators on 

the vacuum states. We derive the boundary S'-matrix and present an integral 
formula for the correlation functions. The boundary magnetization exhibits 
O ' boundary hysteresis. We also discuss the rational limit, the XXX model. 



1 Introduction 



In the standard treatment of quantum integrable systems, one starts with a finite 
cd ' size system and imposes periodic boundary conditions, in order to ensure the 

commutativity of the transfer matrix. Recently, there has been increasing interest 
in exploring other possible boundary conditions compatible with integrability. 

For lattice models with non-periodic boundary conditions, the works on the 
Ising model are perhaps among the earliest (see McCoy- Wu [|l|,^, Bariev [^). 
One should also mention Gaudin's work Q concerning the Bose gas with a delta- 
function interaction and the XXZ model (see also [g]). A systematic approach 
to this problem was initiated by Sklyanin ^ in the framework of the algebraic 
Bethe ansatz. It is well known that, in the periodic case, one can construct a 
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commuting family of transfer matrices from a solution R{C) of the Yang-Baxter 
equation (YBE). Sklyanin showed that a similar construction is possible with the 
aid of a solution i^(C) to the boundary YBE (also known as the reflection equation) 

K2{C2) i?2l(ClC2) i^l(Cl) Rl2{Cl/C2) = R2l{Cl/C2) i^l(Cl) ^12(6(2) i^2(C2) • 

(1.1) 

Several groups have investigated the solutions of (|1.1| ) along with the Bethe ansatz 
equations for the models associated with them [J7|,|^j9|,p^|ri,12,13|. For a recent 



review and further references, see [14|. 

The boundary YBE ( [1.1| ) was originally formulated by Cherednik |15], as a fac- 
torization condition for scattering at a boundary wall. In massive integrable field 
theories, Ghoshal and Zamolodchikov ||l^ discussed integrable massive deforma- 
tions of conformal field theories in the presence of boundaries ll^, and developed 



a bootstrap approach to them (see also [18,1£,2C] for related works). 

One model treated by Sklyanin's method is the XXZ spin chain with a bound- 
ary magnetic field: 



Tj-bare 



1 



00 



2,=i 



(ag+ia^ + al^^al + Aa^+icr^) + haf 



(1.2) 



Here, we consider the model in the limit of the semi-infinite chain, in the ferromag- 
netic regime, A < — 1. In this paper, we present an alternative approach to the 
solution of this model, applying the diagonalization scheme developed in [21,p^] 
for the 'bulk' Hamiltonian 

1 oo 

^bare ^ _ ^ [af^.a^ + ^l^.al + AaU,a{) . (1.3) 



For a systematic account of this method, see [23|. 



The ground state energy of the bare Hamiltonian (1.2) is a divergent scalar, 
c(A, h), which we subtract. We subsequently consider the spectrum of the 'renor- 
malized Hamiltonian,' 

Hb = H%-'^'-c{A,h), (1.4) 



which has a ground state energy equal to 0. 

The starting point in the bulk theory is to identify the space of eigenvectors of 
the Hamiltonian (^]j) ('space of states') with the tensor product 7i 0Ti*, where 
7i = 7i^^' © 7i^^' is the direct sum of level one integrable modules of Uq{sl2), and 
TC* its dual. Here, 



A 



+ 



^-1 



(1.5) 



and -1 < q <0. 

In contrast, we relate the space of states of the boundary Hamiltonian ( |1.2| ) 
to the half space 7i. Accordingly, we rewrite Sklyanin's transfer matrix Tb(C), on 
the semi-infinite lattice, in terms of the type I vertex operators <l'e(C) (see (|2.13| ) 
below). The renormalized Hamiltonian ( |1.4[ ) is defined as the derivative of Tb{C) 

at C = 1- 

It should be emphasized that, unlike in the bulk theory, the algebra Uq{sl2) 



does not commute with the boundary Hamiltonian (1.2), and hence is not a sym- 
metry algebra of this model. Nevertheless, we will show that one can use Ti and 
$e(C) to describe the space of states. A simple bosonization formula is available 
for the spaces Ti^"^' and the vertex operators [p^j2^ . Our basic result here is that, 
in the bosonization language, the vacuum vectors in each sector 7i^^> have the 
remarkably simple form: 

\i)B = e^» \i) , (1.6) 

where \i) is the highest weight vector of 7i^^\ and Fi is quadratic in the bosonic 
operators (see ( |3.4| )). This fact allows us to compute the spectrum of (|1.4[), and 
the spin correlation functions. We subsequently find the following features. 

In the presence of the boundary term, the Hamiltonian ( |1.2D lacks spin-reversal 
symmetry, and the two vacuum vectors, |i)_B, carry different energies. This energy 
difference (see ( |2.28| )) corresponds to the binding energy of the boundary bound 
state of ilF 



Although the expression (1^) makes sense in the bosonic Fock space, there is 
no a priori guarantee that li)^ are well defined for all values of h. We surmise that 
they exist as long as they are regular at g = 0. According to this criterion, \'\)b 
exists if and only \i h < hi or h > hi , and \Q)b exists if and only if /i > —hi 

(2) 

or h < —he ■ Here, the critical fields are 

^^ - -4q ~ 2 ' ""- - -4q - 2 ' ^ ' 

We will show, in Section 3, that our criterion for existence agrees with the singu- 
larity structure of the boundary S'-matrix. 

As in the bulk theory, excited states are created by action of type II vertex 



operators ^^(■?) on \i)B- Using our expression for the vacuum states (|1.6D , we 
deduce that single-particle states obey the relations 

K(0\i)B = MJp^iO %{r') \^)B , (^ = ±) . 

The matrix M^^> (^) is therefore the boundary ^-matrix, analogous to that formu- 



lated in |16| for integrable field theories. 



The bosonic expression for the vacuum vectors enables us to derive integral 
formulas for the spin correlation functions. The boundary magnetization is an 
especially simple case of these correlation functions (see (4.19)). For example, 



s(Okf|0)B 
b{0\0)b 




h 


l-q^l + r 



where 

— 4g 1 — r 

The values of h where both vacuum states \i)B coexist are those where boundary 
hysteresis, or wetting, occurs [||,|,25]. 



The text is organized as follows. In Section 2, we briefly review Sklyanin's 
theory, and formulate the vertex operator approach to the boundary theory. We 
calculate the energy difference between the two vacuum vectors \i)B- In Section 3, 
we find explicit expressions for these vectors, by using the bosonization formulas for 
vertex operators. We then calculate the boundary S'-matrix. Section 4 is devoted 
to finding integral formulas for the spin correlation functions in general, and the 
boundary magnetization in particular. The rational limit, which corresponds to 
the XXX model, is discussed in Section 5. In Section 6 we discuss our results and 
note some open questions. 

2 Eigenvalues of the transfer matrix 

2.1 Sklyanin's formulation 

In this section we recall Sklyanin's results in the context of the XXZ chain on the 
finite lattice, 

-, N-l 

Hb"" = -^ E {^t+i^t + ^l+A + Aal^.al) + h^a^ + h^afj . (2.1) 

k=i 

This will motivate our construction of the transfer matrix for the semi-infinite 
lattice below. 

The transfer matrix is constructed as follows. Below, the reader is referred to 



Appendix A for notation. Let R{C) be the i?-matrix (A.l) of the six-vertex model, 
where C is the multiplicative spectral parameter. We also specify a matrix K^Q, 
corresponding to an interaction at the boundary, which satisfies the boundary YBE 



(|1.1|). Here, we consider only the diagonal solution |[T^ 

Nondiagonal solutions correspond to fields coupled to the other spin components. 
The scalar f{C;r) in ( p^ ) is 



mr) = ^^^:^, ^(,.,) = (^!!:ii^(^!f!i^, (2.3) 

and is chosen so that the i^-niatrix obeys the relations 

K{C)K{C'^) = 1, (Boundary unitarity), (2.4) 

K'ai-q-'C') = E Ra"-l'i-QC^)KS'{0, (Boundary crossing). (2.5) 



These ensure that the transfer matrix ( p. 13 ) satisfies the unitarity and crossing 
relations ( |2l6D , ( ^ ). 



A commuting family of transfer matrices is constructed from R and K as 
follows [|: 

T^''{0=tiVo{K+{C)T{C-')-'K^{C)T{C)). (2.6) 

Here 

T{C) = Roi{C)---Ron{0 G End{Vo'^Vi^---(g)VN) 

denotes the monodromy matrix, Vj are copies of C^, and 

K^{C)=K{C,r.), K+{C)=K{-q-\-';r+Y G End (Fq) , 

where r± are arbitrary parameters. Graphically, the transfer matrix ( |2.6| ) is rep- 
resented in Figure 1^. 

With these definitions, the following statements hold M: 



1. The transfer matrices ( |2.6| ) form a commutative family: 

[r|'^(C),r|°(O] = vcc'. 



2. The Hamiltonian (^Tll) is obtained as 



d 



dC 



Aq 



-^Tr'^TiO ,_ = TtS^^"^ + const. (2.7) 



C=l l-q 



where A is defined by (|1.5|) and 



ft, = ^i^i±^ . (2.8) 

— 4g 1 — r-t 

In Q , an algebraic Bethe ansatz is constructed for the eigenvectors of the transfer 
matrix (|2.6| ). The Bethe ansatz equations for ( ^.Ip were previously derived in [Qj5[ 
using the coordinate Bethe ansatz. 

2.2 The semi-infinite spin chain 

We now consider the XXZ Hamiltonian in the limit of the semi-infinite chain, 
(|1.4| ), with h = h-. and r = r^. Since, under conjugation of Hb by the spin- 
reversal operator n<^f i the sign of the boundary term is reversed, we can restrict 
our discussion to /i > 0, or 

-1 < r < 1. 

The free boundary condition, h = 0, is r = —1, whereas the fixed boundary 
condition, h = oo, is r = 1. When r = (or r = cxd), the Hamiltonian ( |1.2| ) 



formally enjoys C/g(,s/2)^iivariance [26| 



The transfer matrix Tb{C) corresponding to the limit of ( |1.2| ) is depicted in 
Figure ^. It describes a semi-infinite two dimensional lattice, with alternating 
spectral parameters. 

In the naive tensor product space ■ ■ -(SiV 0V 0V, the eigenstates of Hb which 
have finite eigenvalues span a subspace TC, which we call the space of states. Much 
insight about this space can be gained by examining the extreme anisotropic limit, 
g — > 0. In this limit, the Hamiltonian —2qHB scales to 

Hb,o = 2 E(^^+i^fc + 1) + ^(I^^i ' (2.9) 

where a scalar term has been added to ensure that the lowest eigenvalue of Hb,o 
isO. 

When r = — 1, HBfi has two antiferromagnetic ground states, \p^^'), 

Ip'^^^) = ••• «+ f- (8)f+ (8)f-, |p(^^) = ••• (g)f_ (g)i;+ (g)t;_ (git;+. (2.10) 

Accordingly, TC splits into two subsectors, Ti^^' and TC^^', where TC^^' is the span of 
vectors \p) = (Xifcli^p(fc)) called paths, labeled by maps p : Z>i — > {it} satisfying 
the asymptotic boundary condition 

p(A:) = (-!)'=+* for A; >1. (2.11) 



When r 7^ — 1, the eigenvalues of Hba corresponding to \p^^') and \p^^') differ: 

2(1 — r) 2(1 — r) 

Therefore, \p^^') is no longer a ground state. 

The splitting into two sectors TC^^' {i = 0, 1) described above persists for general 
values of r, and — 1 < q < 0. For (7 7^ 0, an eigenvector of Hb is an "infinite 



linear combination" of the paths with one of the boundary conditions (2.11), with 
i = 0, 1. Near r = —1, each sector contains a state |i)_B, unique up to a scalar 
multiple, which has the lowest energy in that sector. We call 1^)^ the vacuum 
vectors, and extend the definition to all values of r by analytic continuation. We 
remark that these analytic continuations are not necessarily physical states for all 
values of the parameters r, q. 

2.3 Vertex operators 

In order to diagonalize ( |1.4| ) for general values of q, with — 1 < g < 0, we follow 
the strategy proposed in [pl[. Figure § suggests (see chapter 4 of ||2^) that we 
identify Ti^'^' with the integrable highest weight module V{A.i) of Uq{sl2), i = 0,1, 
and the horizontal lines with the components $£ *' (C) of the vertex operator of 
type I: 

$(^-''*)(C) : n^'^ -^ n^^''^ V, $(!-*'*) (c) = Y, '^i^'^'^HO «> ve. (2.12) 

E 

In the present case, the algebra Uq{sl2) does not play the role of a symmetry 
algebra, as it does not commute with the Hamiltonian. Accordingly we shall treat 
Tl^^' and <I>e *' (C) simply as objects which enjoy the properties summarized in 
Appendix A, and disregard their representation theoretical meaning. 

The point of using the vertex operators is that they are well defined objects, 
free from the difficulty of divergence. We define the 'renormalized' transfer matrix 

(2.13) 





r«(c) = 




£ 


-'\C')K', 


(o*i;-'' 


Ho , 


(cf. (I2J)). 


Here, 
















^f-'-'HO 


= 'J>i'~ 


"(- 


.-c). 


9 = 




9^)00 



Henceforth we shall suppress the label i for the sectors when there is no fear of 
confusion. The renornialized Haniiltonian Hb is then defined through the formula 
(cf. 0) 

Aq 



d 

dC 



TsiC) ^^= r^.Hs . (2.14) 



Using the commutation relations of type I vertex operators (A. 6), along with 



the invertibility (A. 10), one can show the following properties of the transfer matrix 
dH): 



[TB{C),TBiC')]=0 (VCCO, (2.15) 

Tb{1) = id, Tb{C)Tb{C') = id, (2.16) 

TB{-q-\-^)=TB{C). (2.17) 

These relations reduce, respectively, to the boundary YBE ( |1 . 1] ) , and the boundary 
unitarity and crossing relations (p.4|), ( p.5| ) for the i^-matrix. 

2.4 Energy levels 

Consider now the eigenvalues of the transfer matrix, 

TB{0\v)=t{0\v) , (2.18) 

where \v) is some eigenvector of the transfer matrix. The relations ( |2.15D -( ^.17| ) 
imply 

t(l) = l, t(C)t(r') = l, t{-q-\-')=t{0. (2.19) 

Multiplying <l>e((^~^) to ( p.l8[) from the left, and using the inversion relation 
(A.lOp and the definition of the transfer matrix (|2.13|), we see that the eigenvalue 



problem ( p.l8|) is equivalent to 

j2ki'{c)'^ac)\v) = m Mc-')\v) , (^ = ±) • (2.20) 

e' 

Consider the region of the parameters where 

l<C<l + e, -e<q<0, -Kr <-l + e, (2.21) 

where e > is sufficiently small. Since C > Ij the lowest eigenvalue of the Hamil- 
tonian Hb corresponds to the largest eigenvalue of the transfer matrix Tb{C) (see 
( 2.14| ), and note that q < 0). In view of the study of Hb^o in Section 2.2, we 

assume that 1) the largest eigenvalue of Tj^ {() has multiplicity one, and 2) the 
corresponding eigenvector has an expansion at g = 0, starting from \p^^'). 



Let A(*)(C) denote the eigenvalue corresponding to the vacuum vectors \i)B of 
Section 2.2. Then in the region ( ^I2T|) , A(°)(C) > A^^HC), and A(o)(C) is the largest 
eigenvalue within Ti^^K These statements are not necessarily true throughout the 
region 

1<C, -l<g<0, -l<r<l. (2.22) 

We infer that 

A(°)(C) = 1. (2.23) 



The reason is as follows. In the region (p. 21), this is the largest eigenvalue of the 



transfer matrix. We assume A('')(^) is analytic when 1 < |("^| < \q~^\- Then the 
solution to ( p.l9| ) is unique, and it i s (p.23| ). 

As supporting evidence, we use ( 2. 231) to construct |0)b in Section 3. We then 



use the result to calculate, in Appendix B, the g-expansion of \Q)b (for r = —1), 
up to the order q^ . We find that this agrees with the (/-expansion of the unique 
eigenvector of Hb of the form \p^^') + 0{q), as required in Section 2.3. Therefore 
we conclude that ( |2.23|) is the correct lowest eigenvalue. 
Next we determine A*^"'^^(C). Set 



where Qp{z) is given in ( |A.S| ). We now show that 

A«(C;r) = A(C;r). (2.25) 

We exploit the spin-reversal symmetry to reduce the calculation of |1)b to that 
of |0)b- Let u : 7i^^' -^ 7i^^' be the vector space isomorphism corresponding to 
the Dynkin diagram symmetry (see [p3|). Then 



Noting the relation 



we find that 



$(0,1) (^) ^ = ^im^Q . (2.26) 



a-K{C;r)a- = A{C;r)K{C;r-') , 



V-' r«(C;r) V = A(C;r)Tf (C;r-i) . (2.27) 

If r is close to —1, the eigenvalue 1 of Tg (^;r~^) is the largest. Hence the 



largest eigenvalue of Tg -^(C; r) is given by ( p. 25) ) 



Now, let us compare the two energy levels e^-^'{r) and e'^^(r) of Hb, corre- 
sponding to \0)b and |1)_b, respectively. Parameterize r as 

-(g2)° if -1< r < 0; 
(g2)° if < r < 1, 

where a is real and positive. From ( |2.23| ) and ( |2.14| ), we have clearly 

e(°)(r) = 0. 

Hence the energy difference is simply 

Ae(r) = e(i)(r) - e^°\r) = e^'^\r). 

From (|]2|) and (|1|), we find 

r e(l)sn(2i^'a, A;') if -1 < r < 0; 

^«(-) = ^1^ if < r < 1 (2.28) 

[k'sn{2K'a,k') ^^ ^^ ^ ^ ^ ^' 

where the Jacobi elliptic functions refer to the nome —q = e~'^^ '^ (see 7.5 in [^), 
and 



2K . , ttK' f2K 

t:i"nn nn I / 



e(e) = sinh — ^dn 0, k] , ^ = -ie'" , (2.29) 

vr K \ TT J 

, , 2K]^ ^ ttK' , , 

e 1) = sinh . (2.30) 

vr K 

In Figure |3|, the energy difference Ae(r) is plotted as a function of h. The 
energy increases monotonically from to e(l) (which is the mass gap, see below) 
when < h < he , and from e{l)/k' to oo when he < h < oo. In the latter 
region, the energy of \1)b is greater than the mass gap. The possible existence 
of such a situation, where the boundary bound state energy is greater than the 



single particle energy, was mentioned in |16|, although it does not occur in the 
Ising model. When he < h < he , the state \1)b is well defined, as we shall 
argue in Section 3.1. 

In the next section we shall find the vacuum vectors \0)b and \1)b explicitly 
by using the bosonization formulas for the vertex operators. Once the vacuum 
vectors are found, it is possible to create the excited states by application of the 
vertex operators of type II, ^t(C) (/^ = =t)) i^i much the same way as in the bulk 



theory (see |^^). The commutation relations ( [A.8| ) imply that 

KA^i)---KJU\i)B 

10 



is an eigenstate of Tb(C) with eigenvalue A(*)(C) x IljLi ''"b(CiCj)) where 

rB{C,0=r{C/Orm (2.31) 

and r(^) is given in ( |A.9| ). The single particle energy is therefore given by e'-*'(r) + 
e(^), with its minimum value being e^^'{r) + e(l). Therefore the mass gap is e(l). 

3 Vacuum vectors and the boundary ^-matrix 

3.1 Vacuum vectors 

The perturbative calculation of Appendix B suggests that the vacuum vector \i)B 
is uniquely determined by the relation 

T|\0\^)B = A(^\C)\i)B. (3.1) 

In this section we invoke the bosonization method to find the explicit formulas for 
them, assuming uniqueness. 

First consider the case \0)b- Since the total spin is conserved, it should be a 
(possibly infinite) linear combination of the states created by the oscillators a-n 
(n > 0) over the highest weight vector |0). We make the ansatz that it has the 



form (1.6), with 



n=l L J L J n=l 

Note that the coefficients an, Pn cannot depend on ( due to the commutativity 



of the transfer matrices. Using the bosonization formula ( |A.13 ), equation ( 2.20| ), 
with t{C) = A('^)((^) = 1 and e = —, becomes 

(^(z;r)e-P(^)e<3We^°|0) = { z ^ z''^ ). 

The presence of e^° has the effect of a Bogoliubov transformation 

e-^»a„e^» = a„ + a^a-n + ^^^/3i°\ e-^Oa_„e^« = a_„. (3.2) 

n 

By a straightforward calculation, the coefficients a^ and f3n are found to be 

„7n/2^n „5n/2Q _ „n\ 

an = -Q^^ /?(°) = p^ - en ^^ , 

11 



where 

1 if n is even; 

if n is odd. 

We must verify that the state \0)b found in this manner satisfies also the 

relation (|2.2C| ) with e = +. Using 

the verification is reduced to showing that 

dw {1 - rz){w^ - q^){w - rq^) r-(^)^^-(^6^-i. ^ . -i ^ 

2iTiw{w-q^z){w-q^z)iw-q^z-^) ^' ^ '' 

Here the contour encircles w = Q,q^z but not (f'z . The desired equality can 
be shown by the change of variable w — > q^w~^. 

Similar calculations give the other vacuum vector |1)b, as well as the vacuum 
vectors in the dual space, 7i*^^> , defined by 

B{i\TBiC)=A^'\OB{i\. 

The analog of ( g?20D is 

j2Bmuc')Ki'{c) = A«(c) BiWAo, is' = ±) . (3.3) 

e 

The results are summarized as follows: 

-.CO oo 

H)b = e^^\^), F, = Ij: ^a\ + Y. /?««-«, (3.4) 

n=l I- -I I- -I n=l 

-. oo oo 

b{^ = (^Ie^% G. = \i: r5o«n + E ^^-n, (3.5) 

n=l '^ -"- -I ra=l 



where 



an = -g'^", 7n = -'Z"'", (3.6) 



Injl^n 



[2n] I , g''"/^ r- n 

+ [2n] ^' " ^' 



5n/2^n 



„-3n/2/i _nN j^^-j (^ — 0); 



12 



The relation ( 2.27] ) suggests that the two vacuum vectors should be related by 

u{\0)b) , = \1)b. (3.9) 

However, since the spin-reversal symmetry is obscured in the bosonization, we do 
not know how to verify this directly from the explicit formulas (| 

Remark. It is useful to generalize ( 2.13| ) slightly and consider 



4\C,t) = gY.<^f''-'\tc')Kl' {C)<^l]~'^'\tC), (3.10) 



£,£' 



which corresponds to changing the spectral parameters for the vertical lines from 
1 to t~^. The properties ( p.l5 )-(2.17) are unaffected by this change. The vacuum 



vectors for (|3.1[l|) are obtained by replacing 

an -^ t-'^'^an, a-n -^ t^^'a-n {n > 0) 
in (3.4), (3.5). Doing so and specializing r = (^^ (i = 0) oi r = C,"^ {i = 1) we can 



drop one of the terms in ( |3.1[l| ) , to obtain 

BmUtC^)^-{K) = a(C')b(o|, $*_(tc-i)$„(iC)|o)B = a(c')|o)b, 
B{W+{tC^)^+{tO = a(C')b(i|, ^\{tC^)^+{KmB = A(C')|i)b, 

with 



\{z) 



loo 



iq^z'^;q^)oo (g2z-2.g8^ 



Joe 



These formulas play a role in determining the asymptotic behavior of the iteration 
of the vertex operators, conjectured in [21| and proved recently in p7| . 

Next, consider the small q expansion of the vectors \i)B (i = 0, 1). As stated 
in Section 1, we examine whether they exist or not in the sense they are regular 
at g = 0. Here, we say a vector is regular if it belongs to the upper crystal lattice 
of Kashiwara (see (|B.6| ) ) . 

The expression ( |3.4D for \i)B is not suitable for this purpose, since we do not 
know a simple criterion as to whether a vector expressed in terms of bosons belongs 
to the crystal lattice. In Appendix B, we rewrite \0)b in terms of the global base 
vectors. Although we have no proof to all orders, the expansion suggests that \0)b 
is well defined when /i > 0, but not when \r\ > l/\q\, i.e., —he < h < —he ■ The 
expansion of \'\)b is obtained by the replacement of r to r^^ in ( |B.6| ). Thus, we 
observe three regions for /i > 0: 

(i) < /i< h^}^: Both |0)b and \1)b are well defined; 
(ii) he < h < he : \0)b is well defined, but \1)b is not well defined; 

(2) 

(iii) he < h < oo: Both \0)b and \1)b are again well defined. 



13 



3.2 Boundary S'-matrix 

The two particle 5-niatrix is unchanged from the bulk theory: 

^(0 = -R{CY''' = -RiO- 



The formula ( |2.31| ) shows that the eigenvalues corresponding to the states 
^liiOH) B and ^^('^~^)K)_B are equal. In fact, these states are proportional to 
each other. Using the bosonic expressions for the states b(^| and |i)_B, we find 



KiO\i)B = Ml,^^^{0%{r')\i)B, 



^m^.{c 



B 



m.iOM^^^'^iO- 



(3.11) 
(3.12) 



We therefore interpret M'*' (^) as the boundary S'-matrix. It is computed as in the 

(3.13) 



previous subsection, where we now use ( A. 15 ) and ( A.lq ). We find 
where 



ic2 



M{C;r) = = 



1 



/(e;0 



1 — rq ^ 








(3.14) 



and 



m;r) 



>Mr';r) 



Lp{z;r) 



[q^rz; q 



'^oo («'-'; 9«) 



Notice that, up to a scalar factor, the M-matrix is of the same form as the K- 
matrix ( p. 2D , except for the shift r -^ rq^^. The M-matrix has properties similar 
to those of the if -matrix: 

M2(6)5i2(66)Mi(ei)S2i(ei/6) = ^i2(ei/6)Mi (6)521(66)^2(^2), 

MiOMiC') = 1, 

M'A-Q-'c-') = E s':,~:'i-qe)M^;{o- (3.15) 



These formulas can be obtained either directly, or by using ( |3.11 ), ( |3.12D . For 
instance, the last relation follows if we multiply ( p. 11 ) by ^^i^') from the left and 
take the residue at 6 = C using the property ( A.12| ). 
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In [^], similar equations to ( 3.15 ) were proposed, as the bootstrap conditions 
for boundary field theories. Here, these equations are derived, in the context of 
the lattice model. 

We now consider the analyticity properties of the boundary S'-matrix, whose 
singularity structure determines the existence of the boundary bound states. 



First, recall the discussion of the Ising model (in the continuum limit) in [16|. 
The boundary S'-matrix is given by 

«.(«) = Ua„h(---);-^--j^ (3.16) 

where 6 is the rapidity of the particle with mass m, and k = 1 — h? /Im, h being 
the boundary magnetic field. The result is as follows: If < /i < \/2m, the pole 
of Rh{9) at 6 = ^ — iv, where k = cos v, stays in the physical strip, < Im 6 < tt. 
When h > y/2m, the pole moves away from the physical strip. 

In our case, the pole at ^^ = rq~^ is in the physical strip, 1 < |^^| < \q\^'^, 
when —1 < r < q {0 < h < he ), but also when —q < r < 1 {he < h < oo). 
Here, we considered only the region —l<r<l{h>0). These two regions are 
the regions (i) and (iii) discussed in Section 3.1. 

We recognize that 1 1)^ is the bound state corresponding to the above mentioned 
pole, by using the bosonic expressions for the vacuum and single-particle vectors. 
The single-particle state ^+(C)|0)_b has a series of simple poles at 

f = q^^r; g"V"\ q'^r"'^, ••• ; q^r, q'r, ■■■ . 

By explicit computation, we find that the state \V)b can be obtained as the residue 
at the first pole: 

Res^2=,-i, L(0|*+(r')^ j = 9 c{r-^) b(1|, (3-17) 

Res52=,-i, U;(O|0)b^ j = g c{r-^) |1)b, (3.18) 



where c{z) = {q^z^\q^)oo/{q^z^]q^)oo- The equations ( 3.17[) , (|3.18[) are also valid 






if we replace B(0|*+(r') ^ B(l|*-(r'), **+(0|0)b ^ ^-(0|1)b and 

The formulas above allow us to interpret \i)B as the 'boundary bound states' in 

the regions (i), (iii) mentioned above. 
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4 Correlation functions 

4.1 A^-point function 

In this section we calculate the vacuum expectation values of products of type I 
vertex operators, and obtain them as integrals of meromorphic functions involving 
infinite products. Upon specialization of the spectral parameters, they give multi- 
point correlation functions of the local spin operators of the XXZ chain with a 
boundary interaction. 

We will consider the following A^-point function with N even: 

PkLJCu ...,Cn) = ^^r^^ . (4.1) 

Fixing {ei, . . . , eat}, let us denote by A the index set 

A = {j\l<j<N,ej = +l}. 

Since the total spin is conserved, ( |4.1| ) is non-trivial only if J2j=i ^j — 0) ^^ which 
case A has N/2 elements. 

We remark that, as a consequence of the relation of the two vacuum vectors 
( |3.9| ) and the one for the vertex operators ( 2.26| ), one obtains 



4'.^..e.(Cl, ...,Cn)= P^l...-eJCu •••Xn) __,. (4.2 



In order to evaluate the expectation values ( |4.1| ) we invoke the bosonization 
formulas (|A.13|) , (|A.14D , ( [A. IS]) , ( |A.ig| ) . By normal-ordering the product of vertex 
operators, we have 

N ,,^ / 2 / 4\ 



oo 



W{Wa -Wh){Wa - q 



'w,) 



La-^C, 27r^^ Y[{zj-q-^Wa)Y[{wa-q^Zj) 

j<a a<j 

x/«({z,},{u;4), " " (4.3) 



where each contour Ca is the same one defined in ( [A. 14 ) as Ci. Here, Zj = ("? and 



Tii)(s \ s \\ _ B{i\ew {Y.n=i a-nXn) exp i-Y.'i^=ianYn)\i)B 

b{Ai)b 
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x„ 



7n/2 



N 



7"/2 



E-"-ttE<' yn 



[M ^, ^ [n] ^ 



-5n/2 iV n/2 



\M U ' [' 



n 



Next let us calculate the expectation value l^^'{\zj]^ {wq^). By using the Bo- 
goliubov transformation (p.2|), we obtain 



/^^H{^i},{^a}) 



(z|e^»e-^«exp (^^^^ a_„(X„ - a^i;; 



b(«|«)b 



[2n][n 



exp(-E^/?«^.)exp(iE 



1 ^ f2nlM 



n=l 



2 ^ n 

n=l 



"'^n^n 



We then insert between e ' and e ' the completeness relation of the coherent states 
( C.4 ) and use the integration formula (|C.5| ). As a result we find 



/«({z,},{^4)xn(l-«n7n)^/' 



exp E 



\2n\ \n\ 1 



+ 



fl 2 



^i n 1 - a„7„ 



(4.4) 



b(«|«)b X n(^ -"n7^ 



a/2 



n=l 



= exp lf:MM-^{7./?r + a„<^P + 2/3«5«}). (4.5) 

The sums in the right hand sides are evaluated by making use of the following 
formulas. 



exp f y i!^cz-W2^n^(i) 

n=l 
exp (f: M,-^n^» 



n=l 

oo 



^(Q2z2;g4)^j I [(l-r-iz)(^(z;r-i) 

• l_g2a+5^^2 \V2 J- l_^a+7/2^y^^ 
l_g2a+7^2) I (1 - g'^+3/2r-lu;)-\ 



exp(yi^g3n/2^n^« 



ra=l 



(gV;g^)c 



1/2 



(1 - rz)(p{z;r), 
93(2; r"i), 
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'2n] 



exp > q 

V , n 



.an^n 



5^ 



n=l 



■l_^2a-1^2\V2r ^_^ 



,a-5/2 



rw, 



l_q2a-3^2j I (1 - q^^^/^r-^w)-^ . 



In the above, the upper hne apphes to the case i = and the lower to i = 1. Let 
us summarize below the results of computation. 

The norm of the vacuum vectors: 

b(0|0)b 
b(1|1)b 



(gV;gS)oo 



iq^;q^)oo{q'^r^;q^)oo' 

{q^;q^)oc{q^/r^;q^)oo 



(4.6) 
(4.7) 



The N-point function: 






,6t \ ^ N 1+ 

^] {q'-.q'rj'"^- 

oo / 



nc7 



i-j+Af 



n 



{q^ ZjZk}oo{q^ Zj / Zk}oo{q^ Zkl Zj}oo{q^ / ZjZk}oo 



j^j,{q^ ZjZk}oo{q^ Zj / zk]oo{q^ zkl Zj}oo{q'^ / ZjZk}oo 

ft [r^|]oo[g^^^|]oo[r Alloc [gV^lloo 

xn i du;. ^>^-^"(g''^>g;g')oo(gV^>g;g')o 

iyc«2vrV^ n( 



Z4. 






. «2> 



X 



-q ^Wa) YliWa - 
j<a a<j 

Yliq^'^WaWb] q^)oo{q^Wa/Wb; q^)oo{wh/Wa; q'^)ooiq'^/WaWb; q 
a<b 

YlYiiq'^Zj'^a; q*)oo{q^Zj/Wa] q^)ooiq^Wa/ Zj] q'^)oo{q^ / ZjWa, q'^)oo 



X < 



Hr 



w„ 



,-2 



rwa 



yr {q^rzj;q'^)o^ 
f},{q'rz,-q^)^ , 



(for i = 0) 



1 



r-T- (fori = l) 



(4.8) 



Here, the contour Ca encircles the points q'^^z- (for all j), q'^^Zj {a < j) and 
q^~^^^Zj (a > j), Z = 1,2,3,.. but not the point g^r"^, whereas the contour Ca 
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encircles the point q^r ^ in addition to the same points as Ca does. We have 
also set 

{z}oo = [z] q^, q^)oo, [z]oc = {z; q^, q^)oo 

with {z;p,q)oo = I\j,k>oi^ - zp^q^). 
4.2 Boundary magnetization 



We now specialize the formula ( [4.8D to obtain the correlation functions of local 
operators. In particular, we obtain the boundary magnetization {cf). 

Let L be a linear operator on the n-fold tensor product of the two dimensional 
space V(d- • -(dV. The corresponding local operator C acting on our space of states 
TC can be defined in terms of the type I vertex operators, in exactly the same way 
as in the bulk theory [E3J23] . Explicitly, if L is the spin operator at the n-th site 



0-° = 0-° (g) id (g) • • • (g) id {a = ±,z) 



n-l 



the corresponding local operator £(*)" is given by 

£»± = 4^^(1,...,1), (4.9) 

£W^ = E^;l{l,...,l)-E^U^,...,l), (4.10) 



where 



-E'„'(Cli- ■ ■ ,Cn) 



gu ^ <i>*f'i-)(Ci)-..<er'^(Cn-i)i>:(^-''^)(Cn) 

£i,...,e„_i 

^^'^:;'-'\Cn)^t~-;HCn-i) ■ ■ ■ <i>il-''^)(Ci) (4.11) 



with L = i for even n and l = 1 — i for odd n. 

Thanks to the inversion property ( [A.IOD , the n-point function of the product 
of successively ordered operators of the form ( [4.11 ) becomes particularly simple: 



n X 

The integral formula for this function is obtained from (^]^) by the specialization 

ei , £2 ,...., Gat -^ — El, ....,— En, £n! ■■■ -^^1 
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Cl)C2, ,Cn ^ -q Cl,---,-Q Cn,Cn,- ■ ■ ,Cl- 

In the rest of this subsection we will examine in detail the one-point function, 



M^'\C;r) 



B{ 



)B 



,ii) 



B{i\i)B 



9[P'-^[{-q-\,C)-P?-{-q-\,0)- 
By virtue of the formulas ( ^.20 ), (|3.3| ) and (2^), M^^\Cj1~) has the symmetry 

A^W(C~i;r)=7W»(C;r). (4.12) 

From ( f4.2[) we have in addition 

7W«(C;r) = -7W(°)(C;r-i), (4.13) 



which we shall check directly (see ( [4.15| ) below). 

Specializing the integral formula (|4.8D to N = 2, ^i = —q~^C ^^'^ C2 = C; we 
obtain 

gP^l{-q-\,C) = T9(.z)pt\z,r) 



(4.14) 



where 



(«)/ \ / ^""^ 7-./ \ 

p)J{z,r)= f^^^^ ^^—-^= F{w,z,r), 



F{'w, z,r) 



C^^l 27rV-l 



(l-rz)(l 



wz 



) e,4(g2u;2) 



zil 



^wr) Q„2{wz)Qg2{w/z) 



.(0) 



The contour C__L encircles the points q z, q z with / = 1,2,..., but not the 
point q^r~^ . The contour C]^\_ encircles q z, q z~^ with / = 2,3,..., but not 
q^r^^. (There is no pole at w = q'^z~^.) The contour C^j. encircles the point 
^2^-1 ^^ addition to the same points as C^j. does. With the change of variables 
w — > q^/w we find 



p± (z,r) = -pYiz^^,r^^) + G{z,r), 



(4.15) 
(4.16) 



20 



where 



G{z,r) 



Res^^g2^F{w, z ,r ) 



(1 — rz){l — r/z)- 



QAq'r') 



@g2{rz)eg2{r/z)' 
Using ( ^.161 ), ( |4.12[ ), and noting that G{z,r) = G{z^^,r), we have 
-(A^(°)(C;r) + Al(°)(C;r-i)) 

= g{z) ip\_ {z^ ,r) + p_ {z~ ,r) + p,^ {z,r~ ) + p_ {z,r^ ] 
= 2g{z)G{z,r). 



(4.17) 



The integral for p^. {z, r) can be evaluated explicitly by the calculus of residues, 
making use of the quasi-periodicity property Qp{pz) = —z~^Qp{z). We find 



-M^°\C;r) = g{z)[pf{z,r)+p'^^\z,r)) 



^^ (W(l -rz)(l-r/z) 

1=1 ^ 



rz){l — q^W/ z) 



(4.18) 



The symmetry ( [4. 12]) is now manifest in the expression. 
The boundary magnetization M^'^>{r) in the sector i is 



M^\r) 



i{i\al\i) 



B 



iiAi-) 



B 



7W(*)(l;r). 



Noting (4.13), we have 

-A^W(r) = A^«(r-i) = 1 + 2 V iz^ 

7WW(r)-7W(°)(r) = 2 



2lj,\2 



1=1 



^'-'■,q')oo 



iq'^;q*)lo {q^r-.q^Y^q^/r^q^) 



2^2 

oo 



(4.19) 
(4.20) 



When r takes the values re 

,(1) ^^A f,(2) 



(1) 



(2) 

q and re = q, corresponding to the critical 



fields he and he , the difference ( 4.20 ) is 0, as it should be |2^. At these points, 
the two boundary magnetizations are equal (at the edge of the hysteresis loop). 

Since the ground state is \0)b for h > and \1)b for h < 0, the boundary 
magnetization of the ground state is 



M{r) 



(M^^\r) for |r| < 1, 
\M^^\r) for \r\ > 1. 
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We see that at /i = (r = —1) the spontaneous magnetization A4^^'{—1) 
— A^(^^(— 1) is nonvanishing. Speciahzing ( [4.20|) to r = —1, we find 



A^(o)(_i) = _ (^9')- 



(- 



2. r,2\4: ' 

'OO 



(4.21) 



Note that — A^('')(— 1) is the square of the bulk magnetization p8[ . 

One can check the formula ( [4.18| ) by comparing ( [4.20| ) with the derivative of 
the energy difference Ae(r) = e^^'{r) — e^^'{r) with respect to the field h. One can 



verify directly, by differentiating ( |2.28| ), that the following relation holds: 

(4.22) 



^^"("l = A^(i)(r)-A^W(r). 



dh 
5 Rational limit 

Grisaru et al. [|l^ calculated the boundary 5-matrix for the XXX chain using the 
Bethe ansatz method. To make comparison, let us consider the limit q — > —1. In 
order to get a non-trivial boundary term, it is necessary to scale r near the point 
r = 1. Setting 

r = {q^r 

the Hamiltonian (|l]^) scales to 



-. OO 

rj-bare \ ^ / „x „x Jy Jy i „z „z\ i i,_2 



2.=i 



where h = l/2a. It turns out that we should also scale the parameter ^ in ( 3.14| ) 
as it approaches ^ = \/— T. If we set 



e = V^{q' 



-lA/2 



the formula ( 3.14| ) for M(^; r) can be expressed in terms of the g-gamma function. 
For example. 



Mr(e;r) = (g2)-v^^ 

^q^ [l+ 2 ) Tg-l (i + 2 J ^q' [l " ^2 ) ^l' [^ + 
V [j + 2 ) V [j + 2 ) TgS (^4 + ^— j TgS (^1 - 
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2 

7^Ta 



In the limit g — > — 1, the g-gamma function is the ordinary gamma function, 

and 

M+(A;a) _ V^A + a - i 

Mr(A; a) ~ ^/-iX -a+^' 



2 

Denoting by 7^(A) and /C(A;a) the R- and K-matrices in [|lO| (see eq. (5.17), 
(5.19), (5.30) and (5.36) in ^ 



hmSi2(ei/6) = -alniXi - \2)al, 
hm M(e;r) = -a'IC{X;a). 



Fohowing |29|, p. 99, let us modify the signs and define the limiting states by 



|An, • • • , Al)^„,..,^i;(.)B = n i-l^l) 1™ %A^n) • • • %A^l)\i)B (5.1) 

j:odd 

where ^j = -/^(g^)^'^^^^/^. Notice that 

I - Ai),,;(,)B = lim %A-^^'mB = (-l)^+(i-^')/2 hm %,{Ci')\i)B. 



Taking into account the signs properly, we find that in the basis ( |5.lD our bulk and 

boundary S-matrices become in the limit — 7^(Ai — A2) and — /C(A; a), respectively. 

Thus, in the XXX model, there is only one critical field, which is at a = 1/2, 

and it is the limit of our he ■ The pole in the boundary S'-matrix at A = V — 1(2— a) 



'2 
(2) 

is in the physical strip, < ImA < 1, when h > he ' = 1. The energy difference 
and the boundary magnetization in this limit are: 



lim Ae(r) 



TT 



smTra 



lim7W(°)(r) = -l-2^(-l 



00 2 



1=1 
,2 



{l + aY 



1 + ^(4^)-*'© 



where i^' {x) = {dP /dx^)\ogT{x). 

6 Discussion 

In conclusion, we have obtained the following results for the XXZ Hamiltonian on 
the semi-infinite chain ( |1.2D : (i) the space of states and the renormalized transfer 
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matri x (|2.13| ), (ii) the vacuum vectors ( p.^) and the energy difference between 
them ( |2.28| ), (iii) the boundary S-matrix (|3.14| ), and (iv) the correlation functions 
and the boundary magnetization ([4.19|) in particular. We have also found the 



critical values of the boundary magnetic field (L7) at which one of the metastable 
vacuum states turn to be unphysical. Above the second critical field, one of the 
vacuum states has an energy higher than the mass gap. 

Contrary to the case of the infinite chain, the symmetry algebra corresponding 
to the boundary problem is not Uq{sl2), since the boundary Hamiltonian does not 
commute with it. In this paper, we have bypassed this point to get the results. 
Sklyanin has suggested to us that the relevant algebra may be the one associated 
with the boundary YBE as formulated in his paper ]g]. In order to discuss highest 
weight representations for this algebra, we would probably need a central extension, 



just as was discussed by Reshetikhin and Semenov-Tian-Shansky |3C] in the case 
of the quantum affine algebras. 

Our approach allows us to diagonalize the transfer matrix (|2.6| ) on the semi- 
infinite lattice. However, the corresponding six-vertex model contains alternating 
spectral parameters (see Figure 2). Therefore the corner transfer matrix is not di- 



agonalized by our vertex operators, since in [21,22,23] a uniform spectral parameter 
is assumed. In particular, we have not given a method of diagonalizing the transfer 
matrix acting in the transverse direction. This would involve reformulation of the 
vertex operators to accommodate nonuniform spectral parameters. 

In the field theory limit, there are also two different directions in which we 
can formulate the Hamiltonians, corresponding to the exchange of the space and 
time coordinates. Sklyanin's transfer matrix corresponds to a spatial boundary 
condition. If the roles of space and time are exchanged, the spatial boundary 
condition is replaced by some initial condition in time. In ||l^ this is represented 
by the 'boundary state' \B). Such a state has also been discussed in the context of 



conformal field theory |31,32|. Our expression (|1.6|) for \i)B bears some resemblance 



to the expression for \B), e.g. (3.38) of [16|. However, the two states exist in 



different spaces. The states \i)B correspond (see ( 3.18D ) to the 'boundary bound 
states' of |1£]. 



Although there exists a Bethe ansatz for the model yj5|,^ , the excitation struc- 
ture have not yet been analyzed using this method. In the case of periodic bound- 
ary conditions, only even-number particle states appear in the spectrum of the 
Hamiltonian, although both even and odd states are found using the vertex op- 
erator approach. It is unclear whether all of the states we describe appear in the 
Bethe ansatz for the boundary problem. The boundary scattering matrix ( 3.14| ) 
should also be obtainable using the Bethe ansatz. 

We believe that the case of a non-diagonal i^-matrix, corresponding to a bound- 
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ary magnetic field acting on all three spin components, can be analyzed using a 
similar method to the one described in this paper. Finally, one may ask whether 
there exist difference equations for correlation functions of the boundary Hamilto- 
nian, as for the bulk case. 
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A Vertex operators 

We summarize here some basic formulas for the i2-matrix and the vertex operators, 
following the Appendix in [|3] . 



A.l i?-matrix 








The i?-matrix for the six-vertex model is 






RiC) = 4y 


(1 - CM 

1 - q^-^ 

(1 - ?')C 


(1- 
1- 

(1- 


-q^)C, 

-e)q 



-,2/-1 



e 1 - q'C 



where 



<C) = c 



[qK-hq^)^{qKhq^), 



2a2 



(A.l) 



1/ 



(2;p)oo =11(1 



Z'p 



(A.2) 



n=0 



Let {w+, f_} denote the natural basis oiV = C^. When viewed as an operator on 
V ®V ^ the matrix elements of R{C^ ^"^^ defined by 



^(C) (v, 



V^i 



Y,Ve,<E>Ve,Reli{0- 



£\,£2 
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As usual, when copies Vj of V are involved, Rij{C) acts as R{C) on the i-th. and j-th 
components and as the identity elsewhere. The Yang-Baxter equation satisfied by 



IS 
iil2(Cl/C2)iil3(Cl/C3)i?23(C2/C3) = ^^23 ((2/(3)^13 (Cl/C3)iil2(Cl/C2). (A.3) 



The scalar factor k(C) ( A. 2 ) is so chosen that the unitarity and crossing relations 
are 

i^l2(Cl/C2)i?2l(C2/Cl) = l, (A.4) 

RilliC^/Ci) = RZiii-Q~'Ci/C2). (A.5) 

A. 2 Vertex operators 

We list here the commutation relations for the vertex operators: 

$,,(C2)<l>.,(Cl) = E ReleliCl/C2)<^e'SCl)^e',{C2), (A.6) 

^;;(6)^;^(6) = - E ^w^kei/6)n,(6)^;,(6), (a.7) 

MOKiO = t{c/o%{OMO- (a.8) 

Here 

^(0 = '^"'Tw^i' ®p(^) = (^■,P)oo{pz-';p)ooip;p)oo. (A.9) 

The type I vertex operators satisfy the invertibility 

gY^^UOMC) = id, g<^e,{CWe,{0 = <5.,«id, (A.IO) 

e 

where g is given by 

For the type II operators the corresponding property is 

*t?-^HeiX^-^'^)(6) = Y^^||2(|)* ' +••• (6 -±6), (A.12) 

where 

and • • • means regular terms. 
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A. 3 Bosonization 

For i = 0, 1, consider the bosonic Fock space 

The commutation relations of a„ are 



a.m,a„ = ()m+„,o , [m,ny^O), [n\ = -. 

m q — q~^ 

On the symbols e"^ , the operators e^ , z act as 

g/3_g7 = e/3+7^ ^9.e7 = ^[a,7]e7^ 

where [9, a] = 2, [9, Aq] = and Ai = Aq + a/2. The highest weight vector of W^*) 
is given by \i) = 1 e '. 

We have the following bosonic realization for the vertex operators: 



^r 


"'\) 


*1(^- 


-''Ho 


^f- 


-''Ho 



JCi 2m q(w - q^C^jiw - q^C ) 



f dw (1 - q'^)wC 

tc,2^i q{w - q^e){w - q\^) 

(A.14) 

(A.16) 
X^{z) = e«*(^)e^*(^)0e±-z±^ (A.17) 



where 



oo oo 

O — ra 7r) /9 n ^ / N ^ ^ Or; 



n=l ^^"J n=l ^^^J 

oo oo 

i?±(z) = ±y ^(?T"/v, 5±(z) = Ty ?^9^"/V". (A.19) 



n=l ^ J n=l L ' 

The integration contours encircle w = in such a way that 

Ci : g^C^ is inside and g^C^ is outside, 
C2 : g^C^ is outside and q'^C'^ is inside. 
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These operators act also on the dual (right) modules Ti*^^\ which are defined 
analogously by replacing e '"'""" with e~(^»+"'^). In particular the highest weight 
vector is {i\ = 10 e~ '. The right actions of e^ and z are given by 



eT.e^ = e^+^, e^./ = e^zt'^'^l 



B The g- expansion 

In this section we calculate the ground state vector \0)b of the Hamiltonian (|1.4|), 
as a g-series expansion in terms of paths for r = — 1 |^,^ . Next we develop the 
vector \0)b written in terms of bosonic operators into a similar g-series. We find 
that the two results agree to the order q^^. 

B.l The path expansion 

In the naive approach, we can expand the ground state vector \0)b as a linear 
combination of vectors \p) in the semi-infinite tensor product 0k>i'^k (see ( p.ll )): 



\P) = '^k>lVp(k), 

where 

p:Z>i^{±}, p(A;) = (-l)Mffc»l. 

We call such a map p a path, as well as the vector \p) itself. In what follows, we 
find it convenient to represent a path p by a sequence of non-negative integers 

y={fl,f2r--Jn), fl> f2> ■■■> fn, 

given via the correspondence 

p{k) = p{k + 1) if and only \i k ^Y. 

We also use the convention that |/i, ••-,/«, 0) = |/i, • • • , fn)- The degree of Y is 
/i + • • • + /n- The 'vacuum path' pQ given by po(^) = ("1)*^ is the unique one of 
degree 0, which we represent by the symbol Y = (p. 
Let 

\^)b = \4>) + Y.^'vj (B.l) 

i>i 

be the expansion of the ground state vector for (|1.2D , where e = — 1/(2A) = 
— g/(l + q^). Each Vj is a linear combination of paths 

^3 = E ^A^W)- (B.2) 
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Here we demanded that (j) never appears in Vj (j > 1), by multiplying a scalar to 
\0)b as necessary. 

We are going to solve the following eigenvalue equation order by order in e: 

E li^k+i^l + 1 + ck{e)) - 2eQ | |0)b = 0, 
where 

k>l 

We have included the c-number terms Cfc(e) = J2j>i^'^'^k,j to ensure that the 
eigenvalue is 0. 

Up to the third order, we obtain 

vi = 2\2) + Y,\k + 2,k), (B.3) 

fc>i 

V2 = 6\4)+2j2\k + 4.,k)+2Y,\f^ + 2,k,2)+ ^ \l + 2,l,k + 2,k), 

(B.4) 



fc>l fc>3 k>l 

l>k + 3 



v^ = -4|2) + 20|6) - 3|3, 1) + 2|4, 2) -^\k + 2,k) + b^\k + Q,k) 

k>3 k>l 

+4|4, 3, 1) + 4 E |fc + 4, fe, 2) + 6 E 1^ + 2, A;, 4) + E 1^ + 4, A: + 3, A; + 1, fc) 

fc>3 fc>5 fc>l 

+2 ^ |/ + 2,?,fc + 4,A;) + 2 ^ |Z + 4, /, /c + 2, /c) 

fe>i fc>i 

i>fe+5 ;>fc+3 

+2 ^ |/ + 2,Z,A; + 2,A;,2) + ^ |m + 2,m, / + 2, /, A: + 2, A;). 

(B.5) 



fe>3 fe>l 

!>fe+3 !>fe+3 

m>!+3 



B.2 Monomials in bosons and the upper global bases 



Next, we expand bosonic expression for the vector \0)b ( |3.4D in terms of paths. 
Since \0)b is written in terms of bosons, there is a difficulty in handling its expan- 
sion in the sense of crystal lattice. We expect that each homogeneous component 
of \0)b belongs to the upper crystal lattice (see [^), although we have no proof. 
On the other hand, the global basis vectors have a (/-expansion in terms of paths. 
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(For the meaning of this statement, see |33,21|, in particular, the remark at the 
end of page 99 of [^] and (2^.) In the following we shall relate the vectors created 
by bosons to the global base vectors up to degree 6. This will make it possible (to 
this degree) to develop the homogeneous components of \0)b into a (/-series, which 
can be compared to the results of the previous subsection. 

For a path Y, denote by G{Y) G "H^^) the corresponding upper global base of 
Kashiwara (see Section 2 of |^lj). Up to degree 6, the possible paths and monomials 
in the bosons are as follows. 



degree 


path 


monomials 





</> 


1 


1 


(1) 


e" 


2 


(2) 


O-l 


3 


(3),(2,1) 


a_ie", e— 


4 


(4),(3,1) 


a'^i, a_2 


5 


(5),(4,1),(3,2) 


a„ie-°, al^e", a_2e" 


6 


(6),(5,1),(4,2),(3,2,1) 


a_3, a_ia_2, a^i, e^° 



In the second column we abbreviated X (8) e °'^^ to Xe^ . 

The following is the explicit relation between the two bases: 

pAo 



l<i9e 
.5/2 



G{<p) 
G(l) 

G(2) 

G(2,l) = V 
G(4) = ( 



2c3e'^o+", 



[2] 

,5/2 



pAo 



a_i (XI e 

Ao+a 



Ao-o 



2[2]^ 



r«^i 



^«-2)®e^'», 



G(5) = ( 



a^ + A^«-2We^° 



[4] 



y 2 y 



pAo+o 
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G(4,l) = ^a„i®e'^«-", 



,3 „3 



G(3,2) = (^a^-^a„.)®e^o+« 






g(g2_l) g , </(2</2 + l) 



[6] ^"" ' [4] [2]^-^^-^ ' 6[2]3 -"-i)®« 



G(4, 2) = ,V2 (-^^a_3 + ji^f-^a-i + :i^i:L^±^.3 ^ ^ ^ .Ao 
G(3,2,l) =g6^e^o+2". 



From these data, we can write \0)b (see (3^)) as a linear combination of the 
global basis vectors up to degree 6: 

\0)b = G(0) + grG(2) + gVG(4) + (73^(3, 1) + (7VG(6) 

V[2]rG(4, 2) + AG(5, !) + •••. (B.6) 

The global base vectors are regular and finite at g = 0. We say \0)b is regular to 
mean that each coefficient in the expansion ( [B.6D , regarded as a function of q by 
setting r = ib(g^)", is regular at g = 0. Hence we find, up to this order, that this 
is the case if and only if |r(7| < 1. 
Now, we expand \0)b in q. Let 

G{Y) = \po{Y)) + q\pi{Y)) + q^\p2{Y)) + • • • , 

be the expansion of a global base vector G{Y), where \pj{Y)) is a finite or possibly 
infinite linear combination of paths 

\p,{Y))=J2cj{r,Y)\Y'). 

Y' 

We have, in particular, 

co{Y',Y)=5y',y- (B.7) 

We have the following expansions. 

G{(P) = \(l))+eY,\k + 2,k)+e^(2Y,\k + ^,k)+ J2 \k + 2,k,l + 2,1) 

fc>0 A:>0 l>0 

k>l + 3 
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+e^( J2 \k + 2,k,l + 2,l,m + 2,m) + 2 ^ \l + 2,l,k + 4,k) 



m>0 fc>0 

l>m+3 l>k + 5 

k>l + 3 



+2 Y, \l + ^,l,k + 2,k) + "Y\k + 4,k + 3,k + l,k) + 5Y\k + 6,k) 

fc>0 fc>0 A:>0 

i>k+a 

-J2\k + 2,k))+Oie''), (B.8) 

k>l 

G(2) = (1 - 4^2 + 0(e^)){|2) + e(-|0) + 3|4) + ^ |A; + 2, A;, 2)) 

A:>3 

+s^(9\6) -2|3,1) -^|A; + 2,/e) + 2|4,3,l)+3^|A; + 2,A:,4) 

fe>3 fc>5 

+2^|A: + 4,A;,2)+ ^ |/ + 2, /, A; + 2, fc, 2)) + ©(e^)}, (B.9) 

fc>3 ft>3 

;>fc+3 

G(4) = |4) + e(-3|2) + 5|6) + |4, 2) + |4, 3, 1) + ^ \k + 2, k, 4)) + ©(e^). 



fc>5 

(B.IO) 



Combining ( p.6D -( p.l[l| ), we obtain the g-expansion of \0)b- We have checked 
the equahty between (|B.6| ) for r = — 1 and ( [B.l[) , (|B.3| )- (|B.5| ) up to the third order. 

C Coherent states 

We here summarize formulas concerning coherent states of bosons which are used 
in Section 4. 

The coherent states \^)i and i{^\ in the Fock spaces 7i^^' ,7i*^^' , i = 0,1, are 
defined by 

oo 

10. = exp(X:^-j|-je„a_„)|z), (C.l) 

oo 

.(el = (^|exp(X:pj^Cna.), (C.2) 

where ^n and ^n are complex conjugate parameters. 

Noting that the highest weight states \i) and {i\ are annihilated by the boson 
oscillators a^ and a_„ with n > 1, respectively, one can easily verify 

an\0^=UO^, ^{^W~n = U M (C.3) 

32 



One can also show that the coherent states {l^j} (resp. {i(C|}) form a complete 
basis in the Fock space W^*-* (resp. Ti*^^'). Namely one can verify the completeness 
relation 

,d„,., = / n ^ e-^~" '*I'*-I'|«„>..{&|. (C^4) 

J „ii [n\[2n\ 

Here the integration is taken over the entire complex plane with the measure 
d^d^ = dxdy for ^ = x + iy. In the proof, the following integration formula is used: 

= n (-det A)"'^'exp (1 £ MM^U^^^n), (C.5) 

n=l n=l 

where ^n are invertible constant 2x2 matrices and Bn are constant 2 component 
vectors. 
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{K-)tiC) 



--1 




Figure 1: Transfer matrix on a finite chain. The spectral parameters attached to 
the hnes 1, • • • , A^ are chosen to be 1. 
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Figure 2: Transfer matrix on the semi-infinite chain. By convention, the lattice 
sites are numbered k = 1, 2, 3, • • • from right to left. 
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Figure 3: The energy difference Ae(r) between tlie vacuum states as a function 
of /i, when q = —0.15. The energy a = e(l), and b = e{l)/k' . The critical field 
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